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SUMMARY 
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This report presents information that is required for use of the 
ACCESS- 1 computer program. ACCESS-1 is a research -oriented compu- 
ter program designed to test the actual performance of various new 
concepts and techniques in structural synthesis. The historical 
and technical background of this program is described in detail 
in Refs. 1 and 2, together with a number of well documented numeri- 
cal examples. 

The current version has three types of finite elements, namely, 
truss elements (TRUSS) , isotropic constant strain triangular mem- 
brane elements (CST) , and isotropic symmetric shear panel elements 
(SSP) . Structural weight will be minimized by modifying the sizes 
of these elements — cross sectional areas of TRUSS elements, and 
thicknesses of CST and SSP elements. Design constraints may be 
imposed on nodal displacements, tensile and/or compressive stresses 
in TRUSS elements, von Mises combined stresses in CST and SSP ele- 
ments, together with the minimum and maximum bounds on cross sec- 
tional areas of TRUSS elements and on thicknesses of CST and SSP 
elements . 

There are two distinct general optimization programs which 
can be used in ACCESS-1. One is CONMIN (Ref. 3) which applies a 
modified method of feasible directions and the other is NEWSUMT 
which implements a sequence of unconstrained minimization technique 
(SUMT) , using a modified Newton's method (Ref. 4) for the uncon- 

3 

strained minimizations. 

ACCESS-1 is an all in-core program and all routines are written 
in standard FORTRAN IV language. No advanced coding techniques are 
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used, thus engineers with basic programming experiences can modify 
or restructure the program for their own purposes. Consequently, 
one may make the best use of ACCESS-1 as a research tool to test 
or demonstrate new ideas or techniques through example problems 
of modest size. For example, the basic version of this program 
declares array sizes to accommodate problems having up to 210 ele- 
ments, 70 nodes and 2 load conditions. 
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ACCESS-1 

Approximation Concepts Code for Efficient Structural Synthesis 
Program Documentation and User's Guide 

1 . Introduction 

The ACCESS-1 computer program was developed to demonstrate 
the effectiveness of a structural synthesis capability formed by 
combining finite element analysis techniques and mathematical pro- 
gramming algorithms using an innovative collection of approximation 
concepts. Three types of finite elements are available: namely, 

truss elements of uniform cross sectional areas (TRUSS) , isotropic 
constant strain triangular membrane elements of uniform thicknesses 
(CST) and isotropic symmetric shear panel elements of uniform thick- 
nesses (SSP) . Structures with prescribed configuration and material 
constants are optimized so that their structural weight is mini- 
mized by modifying the sizing of finite elements, i.e., cross sec- 
tional areas of TRUSS elements, thicknesses of CST and SSP elements. 
Design constraints may be imposed on all or selected nodal displace- 
ments , tensile and/or compressive stresses in TRUSS elements , von 
Mises combined stresses in CST and SSP elements, together with the 
minimum and maximum bounds on element sizing variables. 

The fundamental structure of the ACCESS-1 program is illustra- 
ted in Pig. 1. Upon activation, the preprocessor reads input data 
and completes data processing which is not affected by changes of 
design variables (i.e., element sizes), and the results are stored 
in a convenient form for future retrieval. The design process con- 
trol block (DPC) supervises the design iteration procedure as fol- 
lows. The given initial design data is transferred to the approxi- 
mate problem generator, which performs a complete finite element 
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structural analysis, constraint function calculation, as well as 
constraint deletion and sensitivity analysis of retained con- 
straints- Values of retained constraints and their sensitivity 
data together with the closed form description of the weight func- 
tion complete the approximate problem statement which will be 
transferred to the optimization algorithm (OA) block through DPC. 

OA will improve the design using one of the well established con- 
strained function minimization algorithms to operate on the current 
approximate problem statement. An improved design is proposed to 
DPC and the data for structural analysis is updated based on this 
new proposal. This step completes one stage of the design itera- 
tion procedure . The iterative process is then continued until at 
least one of the termination criteria is satisfied. 

Note two important features here. First, the conventional 
finite element method of structural analysis is divided into two 
parts: i.e., the preprocessor and a part of the APG block. It 

is desirable to include as much data processing as possible in 
the preprocessor. Second, the optimization algorithm is asked to 
improve the design with respect to the explicit but approximate 
problem statement, which may or may not be linear. Hence the pro- 
gram used in the OA block may be completely independent of struc- 
tural problems and practically any inequality constrained function 
minimization algorithm can be used. 

Currently there are two distinct general optimization program 
options available in ACCESS-1. One is CONMIN (Ref. 3) which applies 
a modified method of feasible directions and the other is NEWSUMT 


which implements a sequence of unconstrained minimizations technique 
using a modified Newton's method (Ref. 4) for the unconstrained min- 
imizations. Corresponding to each optimizer, a distinct main pro- 
gram and interface subroutines are supplied. The structure and 
program configuration of these two versions are shown in Figs. 2 
and 3. Minor differences between these two versions (especially 
in the DPC block) are attributed to the fact that the NEWSUMT opti- 
mizer includes certain functions associated with the DPC block as 
well as those of the OA block. 

All routines a^e written in standard FORTRAN IV language and 
they have been tested on: (a) the IBM 360/91 using the FORTRAN-H 

and WATFIV compilers at UCLA: and (b) the CDC 6600 using the FTNX 
compiler at UC Berkeley via a remote batch terminal located at the 
NASA Ames Research Center. Implementation on other types of com- 
puters will be straightforward provided those computers have the 
required main memory capacity. The efficiency of ACCESS-1 when 
applied to relatively large scale problems could be improved by 
using advanced coding techniques. However, in its present form, 
it should be possible for anyone with basic practical programming 
experience in FORTRAN to understand and, if necessary, to restructure 
or modify any of the subroutines with relative ease. If new ideas 
or techniques are to be tested, it usually takes a considerable 
amount of time and effort to develop a new computer program. Ex- 
periences show that ACCESS-1 may be used conveniently as the base 
program for the purpose of such experiments. An example of minor 
program modification is given in Appendix E , where a method to re- 
place SSP elements with conventional symmetric pure shear elements 
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is described. 

The ACCESS-2 program, currently being developed, will handle 
significantly larger problems than ACCESS-1, with more involved 
constraints (e.g. thermal effects, fiber composite materials, and 
natural frequency constraints) . The ACCESS-2 program makes effec- 
tive use of dynamic array allocation and auxilliary data storage, 
hence it can solve larger problems than ACCESS-1 using less main 
memory capacity. On the other hand, program modifications of 
ACCESS-2 will require more careful coding and data restructuring. 
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2 . Program Implementation 

Both the CONMIN and NEWSUMT versions of ACCESS-1 may be 
executed as a stand-alone program. All routines are written in 
standard FORTRAN IV language and use only ANS FORTRAN intrinsic 
functions. In additon, if a CPU timing function is available on 
the user's installation, useful CPU time data at the end of each 
stage and also at the end of a job will be printed by replacing 
the dummy routines CPUTIM and CTIME with appropriate ones. Exam- 
ples of these routines are given in Appendix A for the IBM 360/91 
at UCLA, the CDC 6600 at UC Berkeley and the IBM 360/67 at the 
NASA Ames Research Center. 

The declared array sizes of the basic version are determined 
to accommodate problems with 70 elements of each type, 70 nodes, 

2 load conditions and 40 design variables. If desired, the number 
of elements may be increased by using the space allocated to the 
subsequent element types , where the element type sequence is TRUSS , 

CST and SSP. In other words, 210 truss elements may be used, pro- 
vided neither CST nor SSP elements are used, or 140 CST elements 
may be used, if no SSP elements are used. To accommodate problems 
which exceed these capacities, array sizes declared in the program 
must be modified accordingly. 

If program overlay is not used, the basic NEWSUMT and CONMIN 
versions may require as much as 323 K bytes and 366 K bytes of main 
memory on IBM 360/91 at UCLA, respectively. On CDC 6600 computers 
at NASA Langley Research Center, the NEWSUMT version requires 220gK 
words. When program overlay is implemented as given in Figs. 4 and 5, 
the main memory requirements on IBM 360/91 is reduced to 270 K bytes 
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for both versions. Program overlay won't be so effective on CDC 6600 
computers, because the proportion of program instructions in the 
main memory requirement is much smaller than IBM 360 series computers. 

Depending on the problem and/or the choice of parameters, the 
declared capacity of certain arrays may be exceeded; in such cases, 
processing will be terminated automatically, and appropriate 
messages will be printed out. 


3 . Program Organization 

Implementation of the basic procedure outlined in Pig. 1 
is carried out as indicated in Pigs. 2 and 3. Primary functions 
of all subroutines are listed in Table 1, and this facilitates 
understanding of Figs. 2 and 3. Furthermore, all key subroutines 
contain enough comment cards so that the computer program listings 
also serve as a part of the program documentation . 

Data transfer between the subroutines is carried out primarily 
through labeled COMMON blocks. Labeled COMMON blocks appearing 
in each subroutine or function are summarized in Table 2. In case 
it is necessary to modify the array sizes, care should be taken to 
modify all associated array declaration statements. In additon, 
argument lists of the statements to call the following subroutines 
must be modified. 

S ADM 05 SAD007 SAD008 SADMM8 

Some users may wish to improve array allocation efficiency. This 
could be accomplished by allocating a few large arrays dynamically. 
For this purpose the following arrays are suggested as likely candi- 
dates: 

DG: gradient of retained constraints 

AK: master stiffness matrix 

DU: gradient of displacement degrees of freedom 

Note that the selective inverse matrix of AK shares the same memory 
position with DG. 

Currently, two control parameters IDG and IOPT are not util- 
ized for their intended purposes. Therefore, they may be usea to 
provide additional control capability in modifying this program. 


A useful example is given in Appendix D, where creation of an option 
to replace SSP elements with pure shear elements by modifying a 
part of the program is discussed. This example is shown to encourage 
users of this program to modify it, if required, to test new features. 
Another option which is already implemented in the base versions use 
IDG = 5 to remove regionalization of stress constraints. If IDC = 5 
is specified in the input data, stress constraints are imposed on 
all stress-constrained elements, (cee Sec. 4.4) 
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4. Structural Model and Input Data Preparation 

It is assumed that the reader is familiar with elastic structur- 
al analysis via the finite element displacement method, and also with 
associated structural modelling techniques and typical data prepara- 
tion procedure. Sufficient information in preparing the input data 
cards is given in Appendix B, therefore explanations given in this 
section are limited to the subjects which require somewhat detailed 
technical discussion to avoid possible misunderstandings. 

4.1 Node/Element Numbering Scheme 

The solution of linear simultaneous equations is obtained by a 
sequence of calls for SAD007 and SAD008. The coefficient matrix 
(= master stiffness matrix) is stored in a vector form within the 
skyline of the non-zero elements; i.e., there are no operations or 
no storage allocations with elements that remain zero during the 
solution (see Fig. 6 of Ref. 1) . The coefficient matrix is decora- 
posed to LDL form by SAD007 and back anc forward substitutions are 
then carried out by SAD008. The decomposed matrix DL is overwritten 
in the memory area where the stiffness matrix is originally formed. 
The elements of pointer vector IIK indicate the positions of the dia- 
gonal elements for the matrix stored in a vector form. This scheme 
allows somewhat more flexible node/element numbering arrangement 
than the ordinary band equation solver. It is better, however, to 
take the same care in preparing data as for banded matrix solution 
scheme; i.e., differences among node numbers associated with an ele- 
ment must be kept as small as possible for all elements. 

4.2 Symmetric Wing Model 

If the webs of a symmetric wing are modelled with SSP elements, 
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only an upper (or a lower) half of the wing is modelled and x,y 
displacements and loadings are anti-symmetric with respect to the 
x-y plane. Displacements and loadings in the z direction are iden- 
tical for both sides of the x-y plane. For example, if a canti- 
lever beam such as that shown in Fig. 6(a) is to be modelled us- 
ing two SSP elements, then the simplified model should be that 
shown in Fig, 6(b). Note that only half of the loads must be 
applied to the node 3, since the other half is implicitly applied 
to the conjugate node 3' (which does not exist explicitly in 
the model) . The neutral plane coincides with the x-y plane 
and SSP elements are always vertical to the x-y plane. The exam- 
ple 1 given in Appendix C will be helpful in understanding this 
feature. 

4.3 Two and Three Dimensional Structures 

ACCESS-1 treats planar (two dimensional) and spatial (three 
dimensional) structures separately. If a structure is declared 
to be planar by specifying ID = 2, the structure lies on the x-y 
plane and the displacements in the z direction are automatically 
suppressed. In planar structures, nodes whose x and y displacement 
degrees of freedom are free should not be classified as boundary 
nodes . 

4.4 Design Variable Linking 

General concept of design variable linking is discussed in 
Sec. 2.3.1 of Ref. 1. In ACCESS 1, if the sizes of some group of 
finite elements of the same type are controlled by a single design 
variable, these elements are defined to belong to the same linking 
group. Sizes of the elements in a linking group are modified in 
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proportion to the initial sizes given in the input data. 

Also design variable linking groups are used to define "regions 
for the regionalization of stress constraints. General idea of re- 
gionalization is given in Sec. 2.4.1 of Ref. 1. Elements which 
belong to the same design variable linking group form a region and 
only one stress constraint per group and per load condition is 
considered for each group in any stage of the iterative design 
procedure. Selection of the representative elements is not rigid- 
ly fixed/ but dynamically updated at the beginning of each stage, 
if the location of critical stress shifts frequently within a re- 
gion between two consecutive stages, iteration process may be un- 
stable, although this type of instability was not observed in solv- 
ing any of the problems given in Ref. 1. However, if the user 
desires to remove the regionalization of stress constraints, specify 
IDG =5. Otherwise IDG = 0. 

4.5 Configuration/Material Group 

If there are a number of elements of the same type having 
identical configuration and material properties, then these ele- 
ments belong to the same configuration/material group. For exam- 
ple, the single -mate rial planar-truss structure shown in Fig. 7 
has only two configuration/material groups. Configuration/ 
material grouping is used to achieve a reasonable compromise be- 
tween limitations on main memory space and the desire for effi- 
cient run times. The element stiffness matrices in the .local 
coordinate system for unit design variable value are identical 
for all elements in the same configuration/material group. It 
is interesting to note that the local stiffness matrices of CST 
elements are independent of absolute edge lengths and only depen- 







dent on shape. 

4.6 Computation of Constraints and Control Parameters 

All constraints are normalized so that the constraint function 
assumes the values between 0.0 and -1. 0, approximately . 


Stress Constraints 


0.5 (a-a (uj )/(a* u) -a* L) ) < 0 

3 . Si 3 — 

d 


computed stress 

0.5 ( -a) / (a^ U ^ ) < 0 

S cl a. — 

S 


computed displacement 


<*a 


allowable stress 

Displacement Constraints 

i 

fN 


allowable displacer, er. 

S a 


0.5<6-6^ U) )/(6^ U) -6^ L) ) < 0 

D 


sizing variable 

0.5 (6, (L) -6)/(6 -6 ) <0 

d 3 cl “ * 

(U) 


upper limit 

(L) 


lower limit 


Side Constraints 


1.0 - D^D £ 0 
D {L) /D - 1.0 < 0 

As explained in the following section, optimization will be carried 
out in the linked reciprocal variable space. Therefore, for stati- 
cally determinate structures, all constraints shown above are linear 
in this space, including side constraints. 

In the preprocessor (SETCON) , all constraints are identified 
and after deleting strictly redundant side constraints, they are 
enumerated and associated pointer vectors to characterize them are 
prepared. After structural analysis in the APG block, all con- 
straint values are evaluated. Due to constraint regionalization 
and truncation based on the computed constraint values, a significant 
part of constraints are truncated from further consideration during 
the particular design stage. Then sensitivity of these retained 
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small set of constraints are computed with respect to the linked 
reciprocal variables. 

• JSIGNG: sign convention of inequality constraints 

Feasible regions in the design space are defined as follows: 

h {!)) _< 0 q = 1,2,...NTC Structural Analysis 

CONMIN optimizer 

h q (B) > 0 q = 1,2,... NIC NEWSUMT optimizer 

However, the NEWSUMT optimizer has a built-in option to 
accept an analysis program, in which feasible regions are 
defined for non-positive values of design constraints. 

This option is activated by specifying JSIGNG equal to -1. 
•SPM: starting point margin 

If an initial design is infeasible, the initial design is 
uniformly scaled up so that all constraints become satisfied 
with certain margins . The minimum margin for the most criti- 
cal constraint is given by 

Max[h (D) ] = -(SPM-1.0). 

q q 

If it is necessary to change the scaling procedure, the sub- 
routine SUBALY must be modified. 

'TRF,TRFINC and TRFMAX: constraint truncation control parameters 

In the APG block, when all constraint function values hg(D) 
are evaluated, critical and potentially critical constraints 
are selected to form the explicit approximate problem state- 
ment. A constraint h (D) is to be retained as critical or 

SI 

potentially critical if 

h_ (D) > C - TRF* [Max (h (D) ) — C] = TBV 
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where 


Max{h (D) ) : the maximum constraint value in each type of 

constraint 

C-: preassigned constant in SETPOS 

-].0 for stress and displacement constraints 
-1.2 for side constraints 

TRF: initial value is given as an input data and 

modified at the end of each design iteration 
stage by TRF = Min (TRF*TRFMUL , TRFMAX) 

The relation between TRF and the truncation boundary value 
TBV is illustrated in Fig. 8. Note that TBV is gradually 
lowered in absolute value so as to truncate more constraints 
as the design procedure converges. 

4.7 Optimizer Control Parameters - NEWSUMT Version 
•EPSEA, EPSARS, EPSODM: convergence criteria 

EPSEA: Stage convergence criterion 

Iteration process will be terminated, if three consecutive stages 
produce designs which satisfy 

(OBJ p - OBJp^/OBJp < EPSEA 

(OBJ . - OBJ 0 ) /OBJ .I < EPSEA 

p-1 p—2 p-1 

EPSARS: Convergence criterion applied to the results of 3 sequential 

unconstrained minimizations without updating the approximate 
problem statement. This is applicable only if MAXARS > 3. 

EPSODM: Convergence criterion in the golden section minimum 

search procedure. Convergence is achieved if at a certain 
state, 

( |TLL~TL| + |TL-TR| + |TR-TRRj ) / jTR+TL: s. EPSODM 
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and 

ARR - ALL £ 0.05 
are satisfied. (See Fig. 9) 

•STEPMX: Maximum allowable change in any components of the design 

vector during a single stage in the NEWSUMT optimizer. Usually, 
it is not necessary to use this feature, but if some constraints 
are found to be highly nonlinear and errors due to the first 
order Taylor series expansions are excessive, this parameter 
will be useful to confine the design within a reasonable range 
during one stage of the overall iterative design procedure. 

•DELTACs Initial transition point for extended penalty functions. 

4.8 Optimizer Control Parameters - CONMIN Version 

•EPSSTGs Same as EPSEA in S. 4.7. 

•EPSVJK: Same as EPSVJK in §. 4.7. 

* ITMAX : Maximum allowable number of iterations in the CONMIN opti- 

mizer. Here, one iteration is equivalent to one direction 
search followed by a one-dimensional minimization. 

•CTL: Initial width of active region of constraints. A con- 

straint is defined to be 

violated if h g (D) > 0 

active if 0 > h (D) > CTL 

y 

non-active if CTL > h (D) 

y 

Note that CTL < 0. (default value = -0.01) 

•CTLMIN: Upper limit of CTL. This is not an important parameter 

for ACCESS-1, and it is recommended that the default value 
of -0.001 be used. 
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•DELFUN: Convergence criterion among one dimensional minimizations. 

Iteration process will be terminated, if in two consecutive 
iterations, ABS {1. 0 “OBJj_^/OBJj) < DELFUN, and the current 
design is feasible. 

4.9 Printout Control Parameters 

There are two parameters used to control the line printer 
output quantity, namely IPRINT and JPRINT. The greater the inte- 
ger numbers assigned to these parameters, the more detailed out- 
put will be printed. 

* IPRINT controls printouts from all programs except those from 
optimizers. Brief summary of the output items is given in 
Table 3, Standard output will be obtained by assigning 
IPRINT = 2. 

•JPRINT controls output from optimizers (see Table 4) . Standard 
values are 2 for NEWSUNT and 1 for CONMIN. 

4.10 System of Units 

Input data of ACCESS 1 computer program may be prepared in any 

unit systems as long as they are consistent. For example, if the 

units for length and force are decided to be centimeter and Newton, 

2 

respectively, then the unit for pressure or stress must be N/cm . 
Example problems given in Appendix C are presented both in the 
International System (IS) of Units and in the U.S. Customery (US) 
Units. Computer input data examples are prepared using numerical 
values associated with the us Units , simply because all the examples 
were presented originally in various literature in the US Units. 
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5. Restrictions and Limitations 

As explained in the previous sections, the problem size 
which the base version can solve is limited to 
70 elements for each element type 

20 design variables for each type of element, but the total 
should not exceed 40 
70 nodes 

2 load conditions. 

These numbers may be easily modified by changing the sizes of arrays 
declared in the program. However, it is not practical to solve 
large problems using ACCESS-1, even if the computer has large 
main memory capacity. 

The program permits the imposition of upper bounds on element 
sizes as well as lower and upper bounds on positive and negative 
displacements or stresses, respectively. However, this type of 
constraint, when violated, may cause difficulties in convergence. 
This is because these constraints cannot be satisfied by uniformly 
scaling up the design variables. Both optimizers have capabilities 
to start from an infeasible initial design, however, the iteration 
history may be unstable when one or more constraints are violated, 
especially when the NEWSUMT version is used. This shortcoming will 
be eliminated in ACCESS-2. 
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TABLE 1. FUNCTIONS OF ALL ROUTINES 


Main subroutine in the APG block. Organize finite element 
structural analysis, constraint calculation, constraint 
deletion and selective gradient computation. 

Computes objective function and constraint function values 
for the approximate, explicit problem. 

Pseudo load vectors for displacement sensitivity computation 
are assembled by this routine. Given the product of element 
stiffness matrix and corresponding displacement state of the 
element, each component of the product vector is transferred 
to the appropriate position of the global pseudo vector. 

Computes all constraint function values based on the current 
analysis results. 

Compute stress state of the specified finite element. Axial 
stress for TRUSS elements, three plane stress components and 
von Mises combined stress for CST and SSP elements. 

Organizes selective sensitivity computation of retained con- 
straints. 


CPUTIM Dummy subroutine (c.f. Appendix A). 

CST Computes an element stiffness matrix in the global coordinate 

system for a given CST element. 

DIP.CTN* Determines the direction vector in the NEWSUMT optimizer by 
means of modified Newton's algorithm. 

INPUT Reads input data except for the optimizer control cards. 
Writes the input data in a readable format. 

LINK Forms the initial linking table. 

LOAD Assembles load vectors for all load conditions. 

MAIN* Main program. 

MULCST Computes sensitivity of von Mises combined stress of the 
specified CST element, given displacement sensitivity. 

MULSSP Computes sensitivity of von Mises combined stress of the 
specified SSP element, given displacement sensitivity. 

MULTMV Carries out post multiplication of a vector to a symmetric 
square matrix stored in a vector form. Used in assembling 
pseudo load vectors. 

* Subroutines used in the NEWSUMT version only. 
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Table 1. FUNCTION OF ALL ROUTINES (continued) 


OBJVAL Computes structural weight at the beginning of each stage. 

ODM* Carries out one dimensional minimization by means of the 
golden section algorithm. 


PEEOPT 


PRE? 


SAD007 


SAD008 


Check satisfaction of side constraints for the initial de- 
sign. Computes element weight for unit value of sizing 
variable. Computes weight coefficients for the initial 
stage „ 

Identifies representative elements for each linking anu 
configuration group. Computes element shapes, direction 
cosines of the local coordinate axes, and an element stiff- 
ness matrix in the local coordinate system. 

Decomposes real, symmetric positive definite matrix into a 
product of three matrices; i.e., a lower triangular, a 
diagonal and an upper triangular matrices. [K] - [L] [D]{Lj 


Back and forward substitution to solve a system 
[L] [Dj [L] A U = $ for U. 


SADM02 Computes two pointer vectors, JC and I IK. JC indicates 
boundary conditions for the displacement vectors. IIK 
contains the position of the diagonal elements of the mas- 
ter stiffness matrix. 


SADM05 


Assembles the master stiffness matrix, given element stiff- 
ness matrices in the global coordinate system. 

Called by SADMQ4 and performs additions of element stiffness 
matrices in appropriate positions of the master stiffness 
matrix. 


SELGRD 


Same as SADG08 , except for the additional capability to skip 
processing some of the right hand side vectors. 

Identifies displacement degrees of freedom, which may be 
associated at least one of the retained behavior constraints 
(displacement or stress constraints) . 

Assembles the selective gradient vectors for retained set 
of constraints. 


SELINV Computes selective inverse matrix of the master stiffness 
matrix. 


SETCON 


Performs pre-multiplication of the selective inverse matrix 
to the right hand side vectors to obtain selective sensiti- 
vity of the displacement degrees of freedom. 

Identifies all constraints and prepares arrays used in con- 
straint function evaluation procedure. 


Subroutine used in the NEWSUMT vresion only 






Table 1. FUNCTION OF ALL ROUTINES (continued) 


Set up the posture table by deleting constraints which 
are not likely to influence design process at current 
stage . 

Set up the pseudo load vector for the displacement sensi- 
tivity computation. 

Computes an element stiffness matrix in the global coordi 
nate system for a given SSP element. 

Interface between the finite element analysis program and 
the NEWSUMT optimizer. 

Reads NEWSUMT optimizer control parameters and activates 
NEWSUMT optimizer. 

Primary routine of the design process control block (DPC) 
in Fig. 1. Also organizes the NEWSUMT optimizer. 

Picks up components of the displacement sensitivity vec- 
tors to form displacement sensitivity vectors for a given 
element. 

Picks up components of the displacement vectors to form 
displacement vectors for a given element. 

Forms the total function by summing up the objective and 
penalty functions. 

Computes an element stiffness matrix in the global coordi 
nate system for a given truss element. 

Updates the sizing variables of each finite element. 

Print nodal displacement state. 


Subroutines used m the NEWSUMT version only. These routines 
must be replaced by the routines listed in the next page 
to implement the CONMIN version of ACCESS-1. 
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Table 1. FUNCTION OF ALL ROUTINES (continued) 


CONMIN 

Primary subroutine of the CONMIN optimizer. Organization 
of constrained function minimization procedure by means 
of the method of feasible directions. 

CNMN01 

Calculation of gradient information by means of the one 
step forward finite difference. 

CNMN02 

Determination of conjugate direction vector or direction 
pf steepest descent for unconstrained function 
minimization. 

CNMN03 

Solution of one dimensional search in unconstrained 
minimization using 2-point quadratic/ 3-point cubic and 
4-point cubic interpolations, sequentially. 

CNMN04 

Called by CNMN03 and carry out specified interpolations. 

CNMN05 

Direction finding by the modified method of feasible 
directions . 

CNMN06 

Organization of constrained one dimensional minimizations. 

CNMN07 

Called by CNMN06 and carry out specified interpolations. 

CNMN08 

i Special linear programming algorithm with one quadratic 
constraint. 

CNMN09 

Un-scale and re-scale design variables before and after 
evaluation of the objective function. 

MAIN 

Main program of the CONMIN version of ACCESS-1. 

OPTMIZ 

1 

Implementation of all functions of the DPC (Design Process 
Control) block. See fig. 1. 

TAYLOR 

Objective and constraint function evaluations based on 
the current approximate problem statement. 


Subroutines listed on this page are used in the CONMIN 
version only. 
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(sub) routines 
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Table 2 Labeled COMMON Blocks (continued) 


(sub) routines 


BLKA01 

BLKA02 

BLKA03 

BLKA04 

BLKA05 
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• indicates that the associated COMMON blocks must be 
declared in the corresponding routines. 










Labeled COMMON Blocks 
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Table 2 Labeled COMMON Blocks (continued) 


(sub) routines 


CONMIN version 
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indicates that the associated COMMON blocks must be 
declared in the corresponding routine. 


39 







Messages prior to any error termination 
.Final results 

Time and counting statistics of the job 
Input data summary 

Reduced design variables at each stage 
Element sizing variables at each stage 
Weight information at each stage 
Displacement state at each stage 
Scaling-up information (if any) 

Posture table 

Required memory allocation for DG, DU and AK 
Detailed scaling-up data (if any) 

Stress state of all stress constrained elements 

All constraint values 

Weight coefficients 
Gradient of retained constraints 
Updated linking table at each stage 
Basic pointer vectors, JC and IIK 
Load vectors 


Identification of representative elements 
Element data (lengths, direction cosines, 

element stiffness matrices in the local system) 

Master stiffness matrix 

Element stiffness matrices in the global system 
Constraint identification for all constraints 








Table; 4 Optimizer Printout Control - JPRINT 


All messages above the horizontal line corresponding to 
each value of JPRINT are printed. 



j 

NEWSUMT Version j 

— . r * 

JPRINT 

Information printed 

o 

Messages prior to any error termination 

u 

RiiMT control parameters 

1 

Iteration stage summary ' 

2 

Maximum step warning 

Mf 

Direction vectors 


Response surface convergence check data 

J 

Direction search data 

Detailed one dimensional minimization data j 


CONMIN Version 

JPRINT Information printed 


No printing 

Initial and final function information 

First debugging level. Print all of above 
plus control parameters. Print function 
value and X-vector at each iteration. 

Second debugging level. Print all of above 
plus all constraint values, number of active 
or violated constraints, direction vectors 
move parameters and miscellaneous information 
The constraint parameter, BETA, printed under 
this option approaches zero as the optimum 
objective is achieved. 

Complete debugging. Print all of above plus 
gradients of objective function, active or 
violated constraint functions and 
miscellaneous information. 


* Standard value 









APPENDIX A 


CPU TIMING ROUTINES 

1. UCLA IBM 360/91 FORTRAN-H version 

A subroutine CPUTIM is a dummy routine and the function 
CTIME(l) gives the remaining CPU time in seconds. This function 
is not included in the FORTRAN Library, therefore a dataset 
SYS1.CCNFLIB must be concatinated to SYS1.FORTLIB 
SUBROUTINE CPUTIM(T.DT. IP) 

T = 0.0 
DT = 0.0 
RETURN 
END 

2. Berkeley CDC 6600 FTRX compiler 

A subroutine SECOND (T) is in the FORTRAN Library and T is 
CPU time in seconds used by the run . 

SUBROUTINE CPUTIM (T. DT. IT) * 

IF(IT.LE.O) GO TO 100 
CALL SECOND' (Tl) 

T1 = T1 - TO 
DT = T1 - T 
T = Tl 
RETURN 
100 T = 0.0 
DT = 0.0 
CALL SECOND (TO) 

♦Written by Dr. Joseph Mullen, Jr., NASA Ames Research Center 
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RETURN 





r 


END 

FUNCTION CTIME (I) 

DATA T , DT/O . , 0./ 

CALL CPUTIM(T.DT.l) 

CTIME = 1000.0 - T 

RETURN 

END 

3. NASA Ames IBM 360/67 FORTRAN H 

A function INTVAL and a subroutine SETTIM are in the FORTRAN 
Library and INTVAL gives CPU time used since the last call of SETTIM 
in mili-seconds units. 

SUBROUTINE CPUTIM (T. DT . IT) 

IF (IT. LE . ) GO TO 100 
IT1 = INTVAL (0.0) 

Tl = FLOAT (ITD/100.0 
DT = Tl - T 
T = Tl 
RETURN 
100 END 

FUNCTION CTIME (I) 

DATA T ,DT/0 . , 0 ./ 

CALL CPUTIM ( T . DT . 1 ) 

CTIME = 1000.0 - T 

RETURN 

END 
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INPUT DATA DESCRIPTION 

The input data description in the card image format given at 
the end of this Appendix should be referred to in preparing an in- 
put data deck for the ACCESS-1 computer program. Example problems 
given in Appendix C will also be helpful. 

Input Data Cards. 


I. Job description and heading 1(11 , 79A1) 


II 


The first column is used as follows 

0 or blank: ordinary heading card, whose content in 2-80 

columns will be printed in the first part of the 
output list. 

1 : indication of the last heading card. 

2 : request for immediate normal termination of 

this job. 

Any number of cards may be used to describe or to comment 
the job. Note that the last heading card must have "1" punched 
in the first column. Without this, all of your data may be 
regarded as heading cards. 


Job control parameters j ( 315 ) 
IDG 


IPRINT 
I OPT 


not used 

print out control (see §4.9 and Table 3) 
not used 


III. Basic structural data 


(415) 


IN 

ID 


number of nodes 

spatial dimensions (2 or 3) (see §4.3) 
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IBNs number of nodes where boundary conditions are specified. 
INL: number of load conditions. 

(315) 

IDRT(j): number of linked design variable groups for the jth 

element type, (j = 1,2,3) 


j = 1 TRUSS 
j = 2 CST 
j = 3 SSP 

(315) 

XCRT(j): number of configuration/material groups for the jth 

element type (j = 1,2,3) 

(3X5) 


IETP(j): number of jth type elements 


IV. Node coordinates (13, 2X, 3E15.6) 


N : node number 

X(N) : X-coordinate of the node N 

Y(N): Y-coordinate of the node N 

Z(N): Z-coordinate of the node N 


V. 


Element data 

For each element type, the following sequence of cards is re- 
quired. 


(13) 


element type identification 


Note ; The element type identification cards are necessary 


for all element types, even if the corresponding type 
of elements is not used in the structure. 


(215, 3E15.4, 315, 5X, 15) : 


TRUSS and SSP 


or 


(215, 3E15.4, 515) 


CST 




M 


: member number 


IDVRj (M) : design variable linking group number of the Mth 

member in the jth element type 
DVj <M) : initial size of the Mth member 

DVOLj (M) : upper limit of the Mth member size 

DVLLj (M) : lower limit of the Mth member size 

ICVRj (M) : configuration group number of the Mth member 

INODj(l,M): node number of the P local node 

IN0Dj(2,M): node number of the Q local node 

IN0Dj(3,M): node number of the R local node (CST only) 

IVC{M,j) : side constraint code of the Mth element 

-1: lower limit only 

0* non-negativity only 
1: upper limit and non-negativity 

2: both upper and lower limits 

(215)1 : CST only 

IPGR (M , K ) : pressure load group number for the Kth load condition 

on the Mth element; put ”0" or blank if no pressure 
load on any element for the load condition. 

Note ; CST elements require two cards per element, while 

TRUSS and CST elements require only one card per element 


Configuration/Material group data 

For each element type, the following sequence of cards is 
required. 

► 

(13)1 : element type identification 

Note ; The element type identification cards are necessary for 
all element types, even if the corresponding type of 
elements is not used in the structure. 




.... X 


VII. 


VIII 


Note-; If the element type identification is given as a nega- 
tive number, it indicates that the material constants 
for this type are identical for all the configuration/ 
material groups, thus only the material data card for 
the first group is required. 


(6E12.4) 


ASUj (I) 

ASLj (I) 
RHOj (I) 
Ej (I) 
RNUj (I) 


allowable upper stress limit for the jth element 

type, Ith configuration/material group 

allowable lower stress limit 

specific weight 

Young’s modulus 

Poisson's ratio 


Boundary conditions 


(13, 2X, 

315) 



IBD(I) : 

node 

number of the 

Ith boundary node 

IBX(I) : 

Boundary condition 

codes to x,y,z directions 

IBY(I) : 

0: 

free degree of freedom 

IBZ (I) : 

1: 

fixed degree 

of freedom 


. Load conditions 

For each distinct load condition, the following sequence of 
cards is required. 


(215) 


ILLNOD(K): number of nodes subject to non-zero external loads 

for the Kth load condition 
IPLTYP(K) : number of pressure magnitude groups 


(13, 2 X, 3E15.6) 


ILLB(I,K)s node number of the Ith loaded node for the Kth 


48 


load condition 






CLLMX(I,K) : 


Magnitudes of lumped external loads applied to 


CLLMSf (I ,K) j l the Ith loaded node for the Kth load condition, 
CLLMS(I,K):J in x,y,z directions, respectively. 

(E12 .4)| not required if IPLTYP(K) = 0 

CPLM(I ,K) : magnitude of the Ith pressure group for the Kth 

load condition 


Constraints 

(2E15.4) 

SPM: starting point margin^ (see §4.6) 

TRF: initial truncation factor (see §4.6) 


(315) 
ISCT ( J) : 


(1615) : 


ISC (m, J) 


stress constraint code for the Jth element type 

1 : read the stress constraint codes ISC(m,J) 

element by element 

0 : no stress constraint for all elements in the 

Jth element type 

-1: all elements of the Jth type are constrained by 

lower bounds only 

-2: all elements of the Jth type are constrained 

by upper bounds only 

-3: all elements of the Jth type are constrained 

by upper and lower bounds 

required only for the element type whose ISCT(J) is 
positive 

element stress constraint code for the mth element 
of the Jth type 
-1: lower bound only 

0 : 


no stress constraint 




} 

i 



1 : upper bound only 

2 : both upper and lower bounds 

[ (153 I 

1DCT: number of displacement degrees of freedom on which 

finite displacement constraints are imposed 
(315, 2B15.4) 

NA : node number associated with the constrained dis- 

placement degrees of freedom 
JA : direction (x=l, y=2, z=3) 

ID<;(JA,NA): constraint code 

-1: constrained by lower limit only 

0 : no constraint 

1: constrained by upper limit only 

2: constrained by both upper and lower limits 

DISUL(JA,NA) : upper limit 

DISLL (JA,NA) : lower limit 

(15) 

IVCT: side constraint code 

0: no side constraint 

1: apply the code specified on each element description 

card. (IVC(M,j)) 

X. Optimization control 
KtEWSUMT Version 
(7157 

MAXIFS : maximum allowable number of golden section iterations 

MAXODM: maximum allowable number of one dimensional minimiza- 

tions in an unconstrained minimization 
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MAKARS: maximum allowable number of unconstrained minimiza- 

tions for an approximate problem statement 
MAXNAA: maximum allowable number of iteration stages 

JSIGNG: sign convention of constraints 

1 : if feasible region is defined as h (D) > 0 

g = 1, 2 , » . . /Q 


- 1 : 


if feasible region is defined as h (D) <_ 0 

4 

g — 1 1 2 , t • • • tQ 


JPRINTs print out control for the NEWSUMT optimizer (see 
§4.9 and Table 4) 

(8F8.5, F8 .2, E8.IT 


EPSODMs 


EPSARS: 


EPSEA : 

EPSVJK: 
CUTARP : 
STEPMX: 

DELTAC : 


TRFMUL: 
TRFMAX: 
RPMIN : 


convergence criterion among a seguence of one dimen- 
sional minimizations 

convergence criterion among a seguence of unconstrained 
minimizations 

convergence criterion among a seguence of iteration 
stages 

pseudo load vectors truncation criterion 

response factor (or penalty multiplier) reduction ratio 

maximum allowable change of design variable components 

for an approximate problem statement 

initial transition point for the extended penalty 

function 

trunction factor increment ratio 
upper limit of truncation factor 
minimum allowable limit for response factor 
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CONMIN Version 


(15, 4EX5.5) 

MAXSTG: maximum allowable number of iteration stages 

EPSSTG: convergence criterion amont iteration stages 

TRFINC: truncation factor increment ratio 

TRFMAX: upper limit of truncation factor 

EPSVJK : pseudo load vectors truncation criterion 



JPRINT: print out control for the CONMIN optimizer (see §4.9 

and Table 4) 

ITMAX : maximum number of iterations in CONMIN 

CTL : active constraint width (see §4.8) 

CTLMIN: lower limit of CTL (see §4.8) 

CONMIN iteration convergence criterion 


DELFUN 
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(continued from the previous page) 


lUttOTX O hkLTYPO) 


_x . X * X,* 5 

CPLW (- 4 t) 


WMbkW'PCij, X 


XUNQVffl iPLTfP 


mmtW/A 


mMpumzn 


CJPAM(J. 4 2) 


\cplmC 






l_L 

_J OLLMYC 

c\ 1)._ 


! 




' : l 


i ► ! 










i 

1 

-TtOJsf—J . — 


mSQj /|F JCP4T 




X X X £ X 


2i ... 


x . jK_ x x ^ x V £; x 


I ; , Ca ll M2 ( 


or KBQ 


— 

*— 

P XPi\ 

YP(2): 

, i 



l> Coa/v) TJOAf 2 


tf* Xv 


J/ J l 1? *3 11 *3 ^ ^ *f* *7 id i9 jp 51 S3 ->J U » sa s; M ;M0 t) e7 6) 6* 65 tA> 6 ? 6-3 6? ?0 71 7 

(continued to the next page) 


Data cards which must be present in any case. 
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(continued from the previous paqe) 
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t Data cards which must be present in any case . 

* Either one of the two sets of control cards must be present. 

** Data cards for the other jobs may be placed directly after optimizer control cards, 
t Either this card or four blank cards will terminate this run. 
























APPENDIX C 


DATA EXAMPLES 

L.C. 2 



Material : 
Stress limits: 


Aluminum Alloy E = 0.68948 x 10 7 N/cm 2 


P - 0.0027680 Kg/cm 3 


A 2 

Upper limit = 0.68948 x 10'. N/cm 
Lower limit = -0.68948 x 10 q (truss only) 


Displacement limit: 
Side constraints: 


±5.08 cm in Z-direction for all nodes 


Min. area of truss elements = 0.64516 cm 
Min. thickness of CST and SSP = 0.0508 cm 


Load conditions : 


Two distinct load conditions as shown above 


Example 1. 18 Element Wing Box (IS Units) 
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Material ; 

Stress Limits: 

Displacement Limits: 
Side Constraints: 

Load Conditions: 


7 3 

Aluminum E = 10 psi p = 0.1 lb/in 

A 

Upper bounds = 10* psi 

4 

Lower bounds = -10 psi {truss only) 

±2.0 in. in z-direction for all nodes 
Min. area of truss elements 0.1 in* 

Min. thickness of CST and SSP 0.02 in. 

Two distinct load conditions as shown above. 


Example 1 . 


18 Element Wing Box (U.S. Customery Units) 
(Problem 8 in Ref. 1) 



Example 1 18-Element Wing Box 
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Allowable 
Compression , 
Stress (N/cni 
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2 ~ 5 
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Material: 

stress limits : 

Displacement limits: 
Side constraints: 


Aluminum Alloy E = 0.68948 x 10 N/cm 

o = 0.0027680 Kg/cm 3 

Tension = 27579.2 N/cm^ 

Compression = see table above 

0.8890 cm on all nodes and in all directions 

Min. area = 0.064516 cm 2 


Example 2 25-Bar Truss (IS Units) 
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Mode 

Numbers 

• . 

Allowable 
Compression j 
Stress (psi) ' 
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: 1 ' 

35092.0 j 

2 - 5 
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10, 11 i 
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Material; 

Stress. Limits; 

Cross Sectional Area 
Lower Limits: 
"Upper Limits; 

Displacement Limits : 


7 

Aluminum, E=10 
Tens ion=4Q 000 . 0 

0.01 in 2 
none specified 


psi, / =0.1 pci 

psi, (see Table above 

for Compression) 


0.35 in. on all nodes and in all 
directions 


Example 2 25-Bar Truss (U.S. Customer y Units) 
(Problem 5 in Ref. 1) 
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boundary conditions 


load condition 1 


load condition 2 


constraint specification data 


r r r n 


CONMIN optimizer control 
indicates end of this run, 






























Material : Titanium 

Stress Limits: 125,000 psi 

Displacement Limits: Linear envelope with 100.8 inches 

at the tip nodes. 

2 

Load Condition: 144 lbs/ft equivalent. 


Example 3 Delta Wing {u.S. Customary Units) 
(Problem 10C in Ref .1) 
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for the double precision version 
of ACCESS 1 computer program. 
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APPENDIX D 

OUTPUT EXAMPLE — Example 1 

A part of output listing for example problem 1 (18 element 
box wing) is given. The NEWSUMT version was used with printout 
parameters as IPRINT=2 and JPRINT=2 . From the beginning to the 
end of the first stage is listed without deletion. Printouts 
for the intermediate stages (i.e. stages 2 to 8) cure omitted, 
since they are simple repetitions of the first stage output 
as far as the output format is concerned. The final summary 
of the job is included in the end. 
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APPENDIX E 


PROGRAM MODIFICATION TO REPLACE SSP ELEMENTS 
WITH SYMMETRIC PURE SHEAR ELEMENTS 

Element stiffness matrix for a pure shear element can be shown as 

> 


Ik e ] - 



a 

j 

-1 

a 

1 

Et 

i 

-i 

1/a 

-1 

-1/a 

uim 

a 

-1 

a 

1 


1 

1/a 

1 

-1/a 


where 

E: Modulus of elasticity 

v: Poisson's ratio 

t: Thickness of the element 

a: Aspect ratio a/h 

a: length of the element 

b: full depth of the element 

Stress state is 


a x - 0 

a y = 0 

* _ u + u_ - v 

T =, g , p j & - .a 

T '™TV— " Vi— 5 a 


xy ” JTI+vT 


where 


u , v : x, y displacement of p node in local coordinate syst 

r IT 

5^, v^: x, y displacement of g node in local coordinate syst 

Stress constraint is written as 

IVy 1 i T allowable <in P ut data > 

Only modifications to be made are stiffness matrix, stress coraputa 
tion and gradient of stress computation. By tracing the program 




1 

( 


description, it will be obvious that PREP, COMSTR and MULSSP must 
be modified, accordingly. 

Let us decide to use IOPT parameter and assume that all SSP 
elements are replaced by pure shear elements if IOPT = -1. The 
following modifications will be required. 

PREP 

9 lines below statement number 380 

C-E3 (I)/ (12.0E0* {1.0EQ+RNU3 (I) ) ) 

F=2 . 0E0* (1 . QE0+RNU3 (I) ) /AR 

IF(IOPT.EQ.-l) F=0.0 W Insert 

IH=(I-1)*10 

COMSTR 

14 lines below statement number 250 
SXY=(V(1)+V(3) ) *B“ ( V ( 2 ) -V { 4 ) ) *A 
IF (IOPT . EQ . -if GO TO 252 
STRMIS s =SX*SX+3 . OEO*SXY*SXY 
STRMIS=SQRT {STRMIS ) Insert 

GO TO 260 
252 SX=0.0 

STRMIS=ABS (SXY) 

260 CONTINUE 

4 lines from the end of COMMON/BLKA04 
R=1 .0/ (1 . 0+RNU3 (I ) ) 

IF ( IOPTJEQ. - 1 ) GO TO 100 Insert 

ARSXY=A*R*SXY 


TA{6)=3.0*ARSXY 


GO TO 110 


100 C=0.S*R*E3(I) ! 

t 

IP (SXY.LT.0.0) C=-C 
TAn)=B*DCSPX(M) r— — Insert 

TA ( 2 ) =B*DCSPY (M) 

TA{3) =-A 

TA { 4 ) =TA ( 1 ) i 

TA ( 5 ) =TA ( 2 ) 

TA ( 6 ) =A 

110 DSIG=0.0 


It must be confirmed that IOPT is transferred to all of these three 
subroutines: in other words, check if common block BLKA01 is de- 

clared in each of these routines. 

Also the allowable upper stress limits for SSP elements should 


be selected properly. If the distortion energy criterion (von Mise 

1 


combined stress criterion) is used, t .. .. . - = — c . 

allowable yield. 

If conventional shear stress criterion is preferred, T a n owa bi e = 


1 « 

2 yield* 



